Self-consistent model of unipolar transport in organic semiconductor diodes: 
accounting for a realistic density-of-states distribution 
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A self-consistent, mean-field model of charge-carrier injection and unipolar transport in an organic 
semiconductor diode is developed utilizing the effective transport energy concept and taking into 
account a realistic density-of-states distribution as well as the presence of trap states in an organic 
material. The consequences resulting from the model are discussed exemplarily on the basis of an 
indium tin oxide/organic semiconductor/metallic conductor structure. A comparison of the theory 
to experimental data of a unipolar indium tin oxide/poly-3-hexyl-thiophene/Al device is presented. 
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I. INTRODUCTION 

During the last years the problem of an adequate de- 
scription of charge injection and transport in organic 
semiconductors (OSCs) became of great importance in 
view of application of these materials as basic elements 
of electronic devices such as organic field-effect tran- 
sistors, organic light-emitting diodes, or organic photo- 
voltaic cellsi. 

Organic systems are typically wide-bandgap disordered 
semiconductors possessing as a rule a relatively nar- 
row intrinsic density-of-states (DOS) distribution with- 
out sharp band edges typically approximated by a Gaus- 
sian function^. Transport in such systems is generally 
described by models based on hopping of charge carriers 
between localized states, which are disordered in space 
and energy. In addition to the intrinsic DOS distribution, 
disordered OSCs may exhibit deeper localized states orig- 
inating from impurities, or from chemical and structural 
defects. Those states are normally located energetically 
well below the intrinsic DOS distribution and, therefore, 
they are referred to as deep traps. In general, the energy 
distributions of such trap states are often assumed to be 
also of Gaussian shaped, however for some organic mate- 
rials they are assumed to distribute exponentiallj«^i^. 

To simplify the description of transport, a theoretical 
concept is often applied where charge transport is con- 
trolled by carrier jumps to the so-called effective trans- 
port energy level^"®. This characteristic energy is local- 
ized within the intrinsic DOS-distribution and is defined 
as the energy of those sites a charge carrier visits with 
highest probability, independently of its energetic start- 
ing position. This concept allows one to reduce the de- 
scription of hopping transport to a multiple trapping for- 
malisniS. The transport energy plays then the role of the 
mobility edge in disordered materialsAS and may be used 
to qualitatively separate mobile and immobile carriers in 
OSCs. In this concept, charge transport is provided by 
mobile carriers with a constant mobility where all tem- 
perature dependences in the transport description result 
from the Fermi distribution of carriersii^i^ and from the 
temperature dependence of the transport energyi^. An- 



other version of the transport energy approach utilizes 
the effective, average mobility of the whole ensemble of 
carriers. In this case, the mobility reveals, in general, 
dependences on the temperature and on the carrier den- 
sit- 
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By now, the alternative concept of charge transport in 
OSCs has been also developed where all carriers in intrin- 
sic DOS are considered mobile and their drift mobility is 
assumed to be dependent on the temperature, on the 
carrier density and, in general, on the electric field^^^— . 
Nevertheless, as regards the last dependence, it was noted 
in the literature that that the field dependence of the car- 
rier mobility is not strictly necessary, especially at high 
temperaturesiSiSS. The most advanced version of this ap- 
proach, the so-called extended Gaussian disorder model 
(EGDM)^!^"— , has demonstrated good ability to simu- 
late the current-voltage (/- V) characteristics of different 
OSC devices. 

Another important process providing the functionality 
of organic devices is the injection of carriers from elec- 
trodes into the organic layer in view of a low intrinsic 
charge carrier density. However, the description of this 
process remains still controversial. For charge transport 
in bulk OSCs the electric field at the injecting interface 
is often taken equal to zero assuming space-charge limi- 
tation of the currentii"— li^i^i. Numerous papers assume 
finite values of the interfacial electric field and consider 
injection essentially as a single-particle process. Then in- 
jection is simulated as Fowler-Nordheim (FN) tunnelling 
through the surface energetic barrie r^^i^^ or Richardson- 
Schottky (RS) thermionic emission over this barrier un- 
der an applied external field^i^i^i^"— . In fact, the 
field and the charge carrier density at the interface are 
not known because these values depend on the height 
of the injection barrier which in turn is field-dependent. 
Thus, for a proper description of the charge carrier in- 
jection, a self-consistent determination of the field and 
the carrier density at the electrodc/OSC interface is nec- 
essary which would contain both limits of weak (single- 
particle) and strong (many-particle) injection including 
the space-charge limited regime. An attempt of such a 
self-consistent approach was undertaken in Ref?2i where 
drift-diffusion and Poisson equations in the bulk were 
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considered in conjunction with an injection described in 
the spirit of the RS model. However, apphcation of the 
single-particle injection as a boundary condition for the 
many-particle equations makes this approach question- 
able. In view of disordered OSCs, sophisticated exten- 
sions of the RS model with account of possible injec- 
tion into the tail states below the barrier were devel- 
oped^ii^i^'2^1^'2^'^^ which, however, remain essentially 
single-particle injection models. 

Another approach to charge carrier injection includ- 
ing Schottky barrier lowering and space-charge effects 
was recently developed in Refi^ where three-dimensional 
(3D) hopping of charge carriers on sites of a cubic lattice 
with randomly distributed energy levels was considered 
as well as a sophisticated one-dimensional (ID) continu- 
ous model. The site occupancies and the electric field 
were calculated self-consistently by solving the three- 
dimensional master equation and the Poisson equation in 
successive iterations with account of the field-dependent 
injection barriers. Yet, the effect of the individual image 
potential was thereby overestimated because of duplica- 
tion with the mean field deep in the sample. 

Recently, a self-consistent continuous description of in- 
jection in insulating media in terms of carrier densities 
and mean fields was developed^""— based on the match- 
ing of the electric displacement and the electrochemical 
potential at the interface while still accounting for dis- 
creteness of injected charge carriers. The latter aspect is 
important for a wide range of values of injection bar- 
riers and applied voltages where the individual image 
force dominates the injection process. This ID model 
exhibits a plausible crossover from the barrier-dominated 
behaviour at low voltages to the space-charge-dominated 
behaviour at high voltages and reveals a field-induced 
reduction of the injection barrier as well. The barrier 
lowering relates both to the Schottky effect and the volt- 
age drop in the electrodes initiated by substantial in- 
terfacial charge carrier transfer. The model applies di- 
rectly to inorganic crystal insulators and wide-bandgap 
non-degenerate semiconductors as well as to very narrow- 
band insulators and semiconductors as was indicated in 
Ref.'^''. It applies also to wide-bandgap OSCs if the nar- 
row band approximation assuming a negligible width of 
the Gauss DOS is valid. The latter restriction, however, 
seems to fail in many organic semiconductor o^i'^i^^ and 
excludes the possibility of injection into the tail states 
which proved to be essential for disordered semiconduc- 
tors2JL3^^^. 

In the present paper the self-consistent approach of 
Refsj^i^^i^ is extended to account for a realistic DOS 
shape of the organic material. The charge transport in 
semiconductor/OSC/conductor diode structure is mod- 
eled using the transport energy concept in the spirit of 
Refsjii"— and focusing consideration on the infiuence of 
the injection barrier heights and the DOS parameters 
on the /- V characteristics. As an example, a system is 
studied where only holes are injected from an indium 
tin oxide (ITO) electrode into the OSC layer, whereas 



the possible injection of electrons from the metallic con- 
ductor electrode into the OSC as well as subsequent re- 
combination effects are excluded at this stage, for sim- 
plicity. The obtained model will be compared to experi- 
ments of Ref.— where /- F characteristics of poly-3-hexyl- 
thiophene (P3HT) based unipolar diodes were analyzed 
assuming space charge limited current boundary condi- 
tions at the electrodes. 



II. THEORETICAL MODEL 

Let us consider an OSC layer of thickness L sandwiched 
in between a heavily doped semiconductor and a metallic 
conductor electrode. The organic layer is supposed to be 
extended over the space with —L/2 < x < L/2, whereas 
the semiconductor and conductor electrodes are extended 
over the half-spaces with x < —L/2 and x > L/2, respec- 
tively. ITO, being an electron-conducting semiconductor 
with a deep lying conduction band, is considered as the 
hole-injecting electrode^^i^S while Al is considered as col- 
lecting electrode. The band structure of the system under 
consideration is shown schematically in Fig. [T] 

The OSC is characterized by a DOS represented here 
as a superposition of two Gaussian DOSs (see, for ex- 
ample, Refsi^i^iii^). The first, intrinsic DOS represents 
the highest occupied molecular level (HOMO band) and 
the second one describes the spatially and energetically 
distributed trap states (see Fig. [J). The total DOS dis- 
tribution then can be written as: 
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where Et denotes the average trap energy, Oc and Ot the 
widths of the intrinsic and trap DOS parts, P^ and Pt the 




vacuum 
conductor 




FIG. 1: Schematic band diagram of the considered semicon- 
ductor/OSC/conductor structure and the density-of-states 
distribution in its insulator constituent. The assumed hole 
transport under voltage application is also depicted. 
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numbers of intrinsic states and traps, respectively (notice 
that the level E = coincides with the DOS maximum 
of the HOMO band). 

Spatial distributions of the charge-carrier density p and 
of the electric field Fi in the OSC are described by the 
following system of equations: 



dpc (x) 

kT^, — efiiPc (x) Fi [X) 



dx 



dF, jx) 
dx 



P {x) , 



(2) 



(3) 



where j is the position-independent steady-state current 
density, fii denotes the hole mobility in the conductive 
states of the OSC, its static relative permittivity, k the 
Boltzmann constant, e the positive elementary charge, eg 
the dielectric permittivity of vacuum, and T the absolute 
temperature. The total density of charge carriers p (x) — 
Pc (x) -\- pt {x) is the sum of mobile and immobile carrier 
densities, which are defined as 



Pc {x) 



g{E)f^{E)dE 
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Pt (x) = / g [E) fp (E) dE, 

with fp {E) the Fermi-Dirac distribution of holes, 

Ki [x) + e4)i (x) 
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Here, Ec denotes the effective energy of the transport 
level, Ki the electrochemical potential for holes and 0i 
the electrostatic potential in the organic layer. 

Assuming Ki {x)+e(j)i {x)—Ec 3> kT, Boltzmann statis- 
tics for the mobile carriers can be applied and one obtains 



Pc (x) = Nc exp 



-Kj (x) - e(t)i {x) 
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where 



9 {E) exp 



E 
kf 



dE 



yexp 



■ exp 





erfc 




2 {kTf 


erfc 


(Et~Ec 1 


kT 2 {kTf 





^/2kT 



(8) 

can be interpreted as the effective total number of states 
available for mobile carriers in the transport band of the 
OSC, and erfc (z) is the complementary error function-'^. 



Using Eq. ([7]) in combination with Eqs. ([6]) and ([5]), we 

can express the density pt {x) through p^ {x) as follows 



Pt {x) = Pc (x) 



9 (E) dE 



(9) 



Ie^ Ncexp{-E/kT)+pc{x)' 
and, correspondingly, the total carrier density p (x) as 
^ . r g{E)dE 

p[x)^Pc[x)+Pc[x) — / j^,,rr\ I TT- 

Je^ Ncexp{-E/kT) +Pc{x) 

(10) 

The latter equation together with Eqs. 1^ and Q 
presents a system of equations describing the station- 
ary charge transport in the OSC where the appropriate 
boundary conditions have to be applied. 

The energetic differences between the equilibrium val- 
ues of chemical potential in the electrodes far away from 
the electrode/OSC interfaces and the maximum of the 
HOMO band in the organic layer are defined as the in- 
jection barriers for charge carriers (from now on, the 
minus and plus superscripts denote the quantities at the 
interfaces x — —L/2 and x = L/2, respectively). These 
barriers relate to the difference between the electrode 
work functions E^ (see Fig. [T|), 
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Assuming neither surface charge nor dipole layers at the 
electrode/OSC interfaces one can require continuity of 
the electrical displacement and of the electrochemical po- 
tential across the entire device. This continuity results 
in the following nonlinear boundary conditions'^: 



2 



^F 



L 
±2 



a}, 



'it 
(12) 



where 7^ denote the specific conductivities of the 
electrodes, 6 (z) the Heaviside unit step function, 

— [ps {±L/2)]^^^^ the characteristic distance be- 
tween charge carriers near the respective electrode, and 
Xm = {e/W-Keoei [^Fi (±L/2)]}^^^ the distance from 
the respective electrode to the maximum of the single- 
particle Schottky potential barrier—. The Thomas-Fermi 
screening lengths in the electrodes are introduced as 
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with being the static relative permittivities of the elec- 
trodes, p^ and the equilibrium values of the carrier 
density and of chemical potential in the electrodes far 
away from the electrode/OSC interfaces, the latter ones 
calculated with respect to the bottom of the respective 
electrode conduction band (Fig. [Ij. The last term in 
the exponent of Eq. accounts for the discreteness 
of charge carriers and thus determines the range of the 
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injection barriers and field values where the individual 
image forces dominate the injection process resulting in 
the so-called Schottky lowering of injection barriers^, 
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L 



Details about single particle consideration in the mean 
field description provided here can be found in Ref.'^''. 
Notice that the contribution of the current j in Eq. 
can often be neglected since it is very small in all practical 
cases concerning organic semiconductors^^i^^. 

The nonlinear differential equations ^ and ^ with 
p {x) from Eq. (jlOp and with the boundary conditions 
dm have to be solved numerically. Knowledge about 
the spatial distribution of the electric field gives access to 
the voltage drop V across the system for a given current 
density j, which follows by direct integration of the field 
over the device thickness^i^ and reads: 



V = l+p 



+ 1 



TF 
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L/2 



F, (x) dx - Vbi, 
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(15) 



where —Vbi is the voltage drop in the case of j ~ 0, 
i.e., the built-in potential, given by the difference of the 
electrode's work functions, 
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(16) 



Now, /- F characteristics of the structure under consider- 
ation can be calculated. 



III. COMPARISON WITH EXPERIMENT AND 
DISCUSSION 

To test the presented model, it is applied to experimen- 
tal /- V characteristic&i^ measured on a unipolar device 
consisting of a single P3HT layer of thickness L = 125 nm 
sandwiched between ITO and Al electrodes. It was de- 
duced from different experimentsi^ii^'^ that the total 
DOS in P3HT consists of the superposition of two Gaus- 
sian peaks, similarly to Eq. ([T]). 

Figure [2] depicts I-V characteristics of the 
IT0/P3HT/A1 structure from YleiM measured at 
different temperatures as well as best fits calculated with 
the OSC model parameters indicated. All parameters of 
the electrodes as well as the typical value of the OSC 
relative permittivity, = 3, are taken over from Ref4i. 
At first, fitting the T — 300 K curve, we establish the 
best values for DOS parameters and carrier mobility. 
Next, by fitting of the curves for the other temperatures, 
these model parameters remain unchanged and only the 
barrier heights and the transport energy Ec are 
allowed to vary with temperature. One can see that the 
calculations reproduce satisfactorily the magnitude and 
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FIG. 2: (Color online) Measured I-Y characteristics (solid 
curves) of IT0/P3HT/A1 diode at diflferent temperatures 
and their fitting (dashed curves) with the OSC model pa- 
rameters: [li = 6.7 X lO"* cmV(V s), Pc = 10^^ cm"^, 
Pt = 5x10^^ cm~^ (Tc = 0.035 eV, crt = 0.02 eV, Et = 0.5 eV. 



the general form of the /- Y characteristics in the wide 
range of applied voltages and at temperatures between 
160 and 300 K with the fitting parameters indicated in 
the figure. Remarkable deviation of the theory from the 
experiment in the low- voltage parts of /- F characteristics 
(below —Vhi) is explained by a parallel leakage current 
which is not taken into account in the presented model. 

The obtained fitting values of the transport energy are 
in a qualitative agreement with the temperature depen- 
dence of this value predicted in the literature (see, for 
example, RefSfiii^). However, this energy is not the only 
parameter needed to be varied with temperature. The 
best fitting can be achieved only if we assume the bar- 
riers to change with temperature too. There are at 
least two reasons which could make such changes likely. 
One can see that the ITO/OSC barrier increases from 
0.1 to 0.3 eV with decrease of the temperature from 
300 to 160 K. Such a change corresponds qualitatively 
to the temperature dependence of the bandgap energy 
well-known in classical semiconductors'*^ and recently es- 
tablished in charge transfer complexes based on OSC 
thin filmsi^i^ (it should be noted, however, that such 
an effect is not yet known in P3HT). On the other hand, 
the direction of the barrier change for the collecting elec- 
trode is different from that for the injecting one. It first 
slightly decreases and afterwards does not change with 
the temperature decrease, contrary to the increase of the 
injecting barrier. It may be supposed that the barrier 
heights are additionally affected by dipole layers which 
may emerge at the electrode/OSC interfaces, typically 
for these system o^^'^^ . They can have an intrinsic, indi- 
vidual temperature behavior and thus may contribute to 
the temperature dependence of injection barriers. 

It should be noted that the obtained fitting clearly em- 
phasizes the essential role of the self-consistent bound- 
ary conditions in the proposed theoretical model. In 
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Refi^, experimental I-V curves have been simulated 
within the basically similar mobility edge concept but us- 
ing the temperature-independent transport energy Ec = 
and the space-charge limited boundary conditions of 
Fi (±L/2) = 0. There, the simulation has revealed less 
satisfactory agreement between the calculated and exper- 
imental curves and unreasonably low number of intrinsic 
states of the order of 10^* cm~^ in contrast to the number 
of 10^^ cm"'^ obtained in our simulations. 

In conclusion, a one-dimensional mean-field model 
has been presented which describes self-consistently 
the charge-carrier transport across a semiconduc- 
tor/OSC/conductor structure accounting for effects of 
discreteness of injected carriers and a realistic DOS dis- 
tribution in the organic layer. The energy distribution 
of DOS was modelled by the superposition of two Gaus- 
sian peaks for the HOMO-level and deep traps. In the 
framework of the transport energy concept, assuming 



temperature-dependent injection barriers, good qualita- 
tive and quantitative agreement between the experimen- 
tal and simulated /- V characteristics has been obtained 
in a wide temperature range. Compared with the popular 
EGDM model our approach has a comparable number of 
fitting parameters and gives similar fitting accuracy (cf., 
for example, with Refs.— 
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